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ABSTRACT 



Aims. We present time dependent chemical models for a dense and warm O-rich gas exposed to a strong far ultraviolet field aiming at 
exploring the formation of simple organic molecules in the inner regions of protoplanetary disks around T Tauri stars. 
Methods. An up-to-date chemical network is used to compute the evolution of molecular abundances. Reactions of H2 with small 
organic radicals such as Ci and C2H, which are not included in current astrochemical databases, overcome their moderate activation 
energies at warm temperatures and become very important for the gas phase synthesis of C-bearing molecules. 
Results. The photodissociation of CO and release of C triggers the formation of simple organic species such as C1H2, HCN, and 
CH 4 . In timescales between 1 and 10 4 years, depending on the density and FUV field, a steady state is reached in the model in which 
molecules are continuously photodissociated but also formed, mainly through gas phase chemical reactions involving H2. 
Conclusions. The application of the model to the upper layers of inner protoplanetary disks predicts large gas phase abundances of 
C2H2 and HCN. The implied vertical column densities are as large as several 10 16 cm" 2 in the very inner disk (< 1 AU), in good 
agreement with the recent infrared observations of warm C2H2 and HCN in the inner regions of IRS 46 and GV Tau disks. We also 
compare our results with previous chemical models studying the photoprocessing in the outer disk regions, and find that the gas phase 
chemical composition in the upper layers of the inner terrestrial zone (a few AU) is predicted to be substantially different from that in 
the upper layers of the outer disk (> 50 AU). 
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1. Introduction 

Protoplanetary disks (PPDs) represent an intermediate stage in 
the evolution from dark clouds towards planetary systems. A de- 
tailed characterization of the physical and chemical conditions 
in such disks is of great interest a s they provide the initial con- 
dition s for planet formation (see Naiita et al. 2007; Bergi n~et~ai] 
|2007| for recent reviews on the subject). Information about the 
disk chemical composition has mainly come from radio obser- 
vations, which are s ensitive to the outer cool disk , from ~ 100 
AU up to -900 AU dDutrev et all fl997t iThi et all 12001 . Only 
recently the inner regions of PPDs have been probed by means 
of infrared observations revealing the pres ence of warm C2H2 
and HCN with la rge gas phase abundances dLahuis et all [2006; 
Gibb et all 12.0071) and of complex orga nic species such as PAHs 
dGeers et all 120061 lHabart et alll2006l) . 

Chemical models have mostly concentrated in the study of 
the outer disk regions (e.g. lAikawa & Herbstll 19991) . in part mo- 
tivated by radio observations of several molecular species. Only 
a fe w models have focus ed on the inner disk (< 10 AU) chem- 
istry. IwnJacxetal] d 1 998b studied the chemistry in the disk mid- 
plane, shielded from stellar and interstellar UV photons, and 
found that some organic species could be easily formed on grain 
surfaces jMarkwick et all d2002l) studied the chemistry in the in- 
ner 10 AU including stellar and interstellar UV radiation, stel- 
lar X-rays and radionuclides decay as sources of ionization, and 
adsorption/desorption on grains but not mantle chemistry, i.e. 
formation and destruction of molecules occurs only in the gas 
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phase. They calculated large abundances of organic species such 
as C H4 but they did no t discuss the main chemical r outes to form 
them. ISemenov et all d2004l) and lllgner & Nelson! d2006l) exten- 
sively discussed the ionization d egree over a wide r a nge of radii 
and heights over the midplane. Woods & Willacy (2007) pre- 
dicted that benzene formation can be efficient in the midplane 
of the inner disk (< 3 AU), a dense region protected against UV 
photons where grain surface reactions play an important role in 
building-up a complex organic chemistry. 

In this Paper we investigate gas phase routes to form simple 
organic molecules, such as acetylene (C2H2), hydrogen cyanide 
(HCN), and methane (CH4), in a dense and warm oxygen-rich 
gas exposed to a strong far-UV (FUV) field (hv < 13.6 eV). We 
apply this model to the inner region (< 10 AU) of a protoplane- 
tary disk around a T Tauri star and compare with the abundances 
recently derived from observations. 



2. Chemical model 

Physical models of protoplanetary disks suggests a flared-up 
structure where both the gas density and temperature vary enor- 
mously depending on the radiu s r from the star and heig ht z over 
the disk equatorial plane (see iDullemond et al J 120071 for a re- 
view). The gas density decreases radially outward as a power 
law of r, and in the vertical direction it decreases as z increases 
in roughly an exponential way. The gas heating across the disk 
is dominated by the stellar irradiation of the surface layers at 
large radii (r > 10 AU), thus the kinetic temperature increases 
with z, while at small radii (r ~ a few AU) viscous dissipa- 
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Table 1. Species included in the model 
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tion may beco me an important heating me chanism in the mid- 
plane regions dD'Alessio et aUfl998l fl999). A flared disk is sig- 
nificantly exposed to the strong UV field from the central star. 
The extinction of stellar UV radiation is large in the midplane, 
but it rapidly decreases as z increases. Therefore, there exists a 
layer of low Ay where the disk material is being photoprocessed 
dAikawa & Herbstl 1 19991: IWillacv et all 12000 1. Here we focus 
on the chemistry that takes place in the photon-dominated region 
(PDR) of the inner disk (r < 10 AU), where typically tempera- 
tures are several hundreds of K and gas densities range from 10 6 
to 10 11 cirT 3 . 

In order to qualitatively understand the chemical routes to 
form simple organic molecules in the PDR of the inner disk, we 
have firstly performed a few time dependent chemical models 
in which chemical abundances evolve under fixed physical con- 
ditions representative of this region. We consider an O-rich gas 
with all the hydrogen initially as H2, all the carbon as CO, the 
oxygen in exces s as H2O and all the n itrogen as N2. We use so- 
lar abundances (lAsplund et al.L 120051) and assume that 50 % of 
C and 65 % of O are in carbon and silicate grains respectively. 
Therefore, the initial abundances relative to the total number of 
H nuclei «h, where «h = 2 n(Yli) + n(H), are x(CO) = 1.25 X 
10~ 4 , jc(H 2 0) = 3.5 x 10~ 5 , x(N 2 ) = 3.0 x 10~ 5 , and x(He) = 
8.5 x 10~ 2 . We adopt a gas density of «h = 2 x 10 s cirT 3 and 
run several models with five different kinetic temperatures: 7* = 
1 00, 300, 500, 750 and 1 000 K. 

We first consider a FUV illuminated model with a FUV field 
strength o f y — 50,000 (relative to the Draine interstellar radia- 
tion field; Draine1 [T978l) . the value at 10 AU according to FUV 
observations of various T Tauri disks dBergin et alll2003l l2004). 
For this FUV field and gas density, photoprocessing of the gas 
occurs in a range of visual extinctions Ay ~ 0.1 - 5. If Ay < 
0. 1 the gas is very exposed to the UV radiation and molecules 
are destroyed, and in the case of Ay > 5 photoprocessing occurs 
only marginally. Here we investigate the chemical evolution at a 
mean extinction value of Ay = 2.5. We then consider a separate 
X-ray illuminated model, in which the gas is solely exposed to 
X-rays but not FUV photons, aiming at evaluating the influence 
of these two energy sources on the chemistry. There is evid ence 
for X-ray emission in PPDs (Feig elson & Montmerlel[T999h . Its 
main effect on the chemistry is the ionization of the gas produc- 
ing high energy photoelectrons (ph-e~) which further ionize and 
dissociate the species. To simulate this effect we have enhanced 
the cosmic rays ionization rate by a factor of 1000 over the stan- 
dard interstellar value, thus we take £ = 1.2 x10 s -1 , which 
is rea sonable for the inner regions of PPDs dlgea & Glassgold. 
119991) . 

Th e chemical code used has been described in ICernicharol 
d2004l) . The species included in the model are given in Table [T] 
They were selected on the basis of previous calculations includ- 
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Fig. 1. Evolution of C2H2, HCN and CH4 abundances for a 
chemistry driven by FUV photons (left) and by X-rays (right). 
Different curves correspond to different gas temperatures. 



ing a larger number of species which allowed us to identify the 
most important chemical processes forming C-bearing species . 
The chemical network used was taken from Cernicharo (2004) 
and|Agundez & Cernicharol (l2006). and was checked against the 
NIST Chemical Kinetics DatabasqH an d the latest version o f the 
UMIST database for astrochemistnH (IWoodall et aUl2007l) . We 
include neutral-neutral and ion-molecule bimolecular reactions, 
three body processes and thermal dissociations. 

Important reactions involved in the built-up of small organic 
molecules at high temperatures are those of H2 with radicals 
such as C2 and C2H These reactions have been studied in the 
laboratory over a relatively wide te mperature range: C2 + H2 — > 
C 2 H + H between 295 and 493 K dPitts et allfl982h and CoH + 
H? -> C2H2 + H in the range 178-440 K dPeeters et all Il996t 
Opa nskv & Leonel Il996l) . They have moderate activation barri- 
ers, of about 1400 K, which make them very slow in the cold 
interstellar medium, although at high temperatures they become 
rapid enough to control the abundance of C-bearing species. 
Their rate constant expressions are not included in current as- 
trochemical databases, which are somewhat biased t oward low 
tempe rature chemistry, but can be found in Table 1 o flCernicharol 
(2004). The formation of H 2 on grain surfaces is included with a 
rate constant of 3 x 10~ 17 cm 3 s _1 . 



In order to properly take into account FUV line self- 
shielding, we explicitly calculated a grid of H 2 and CO pho- 
todissociation rates and C photoionization rates, appropriate 
for the physical conditions in the inner disk. This calculation 
was p erformed with the Meudon PDR code dLe Bourlot et all 
Il993l). The model has been described in detail e lsewhere 
dLe Petit et ail 120061: iGoicoechea & Le Bourlol 120071) . In par- 
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Fig. 2. Scheme with the main synthetic routes for the formation of C2H2, HCN and CH4 from C, C + and N. Reactions with a high 
activation energy (Ae) are indicated by a thick arrow. 



ticular, we computed the depth dependent FUV radiation field 
for the range of parameters (x, 7* and initial abundances) 
considered in the FUV illuminated model. H2, CO and C 
photorates were then consistently integrated over these fields 
and used afterwards in the time-dependent chemical calcu- 
lations. D ust properties are known to infl uence the resulting 
rates (e.g. iGoicoechea & Le Bourlolll2007t) . Here we have as- 
sumed the standard g rain size distribution proposed for the ISM 
dMathis et all 1 19771) All the other photorates included in the 
FUV illuminated model have been taken from the UMIST 2006 
database. In the X-ray illuminated model we have included reac- 
tions induced by X-rays, for which we have adopted the cosmic 
ray induced reactions rates enhanced by a factor of 1000. 

3. Chemical routes to C 2 H 2 , HCN and CH 4 

Left panels in Fig.Q]show the evolution of C2H2, HCN, and CH4 
abundances for several gas temperatures in the FUV model. The 
synthesis of these organic species can be divided into three steps. 

(i) Photodissociation of CO and N2 with release of atomic 
C, C + (produced by further photoionization of C) and atomic 
N. The time scale of this step depends only on the FUV field 
strength, thus on^ and Ay. 

(ii) Atom— >molecule transition, from the primary species C, 
C + and N to simple C-containing molecules. Important reactions 
which drive this transition are (see Fig. [2j radiative associations 
such as C + H 2 -> CH 2 + hv and C + + H 2 — > CH+ + hv, and rapid 
neutral-neutral reactions such as C + NO — * CN + O, where NO 
comes from the reaction between N and OH. These reactions are 
nearly temperature independent, thus their rate depends basically 
on the gas density: the higher the gas density the faster they pro- 
ceed. Bimolecular reactions such as C + + H2 — > CH + + H and C 
+ H2 — > CH + H begin to overcome their activation barriers at 
temperatures above 400 and 700 K respectively. Three body re- 
actions are not competitive, compared to radiative associations, 
for gas densities below ~ 10 13 cirT 3 . 

(iii) Processing of simple molecules, CH2, CHJ and CN, into 
more complex species, represented by C2H2, HCN and CH4 in 
our model. As Fig. [2] shows this processing occurs through bi- 
molecular reactions involving H2, C and N. Reactions of neutral 



species with H2 have energy barriers higher than 1000 K, which 
introduces a strong temperature dependence. This is very marked 
for CH4 because its formation involves the reactions with the 
highest energy barriers, making it abundant only at temperatures 
above ~ 700 K. The synthesis of organic molecules in the O-rich 
gas requires that atomic carbon, produced by the dissociation of 
CO, incorporates into C-bearing species faster than reverting to 
CO. This is achieved by different mechanisms at low and high 
temperatures. Below ~ 400 K atomic oxygen is not converted 
into OH (the reaction O + H2 — > OH + H has an activation en- 
ergy of ~ 5000 K) and the main CO-forming reaction, OH + 
C — > CO + H, is inhibited. Thus atomic carbon can react with 
other species instead of reverting to CO. The main obstacle to 
form C-bearing species at low temperature is the energy bar- 
rier of several reactions. For example, C2H2 and HCN reach low 
abundances at 100 K because the reactions of C2H and CN with 
H2 have activation barriers. Above ~ 400 K atomic oxygen is ef- 
ficiently converted into OH, which may react with C to form CO 
but reacts faster with H2 to form water. Thus, most of the oxy- 
gen forms H2O, and CO does not reach its maximum abundance 
allowing atomic carbon to form C-bearing molecules. At high 
temperatures most of the reactions shown in Fig. [2] overcome 
their energy barriers producing a rich C-based chemistry. 

The chemistry driven by X-rays (see right panels in Fig. [TJ 
is rather similar to that initiated by FUV photons except that the 
dissociation of CO and N2, step i, is carried out by collisions 
with photoelectrons: CO + ph-<?~ — > C + O and N2 + ph-e~ 
-> N + N and by reactions with He + : CO + He + -> C + + O 
+ He and N 2 + He + -> N + N + + He. For our chosen FUV 
and X-ray field strengths, these reactions are about 100 times 
slower than photodissociations, which makes the chemical time 
scale to be increased by roughly the same factor, i.e. molecules 
form later (see Fig . |T} . On the other hand, the lower dissociating 
strength of X-rays, compared to that of the FUV field, results 
in a higher steady state abundance of C-bearing molecules. This 
occurs because these steady state abundances depend on the bal- 
ance between the rates at which molecules are formed (through 
chemical reactions) and are destroyed (by either FUV photons 
or X-rays). The destruction rates are lower in the X-ray model, 
compared to the FUV model, and thus chemical reactions have 
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more time to form molecules until steady state is reached. As 
a consequence, in the X-ray model C2H2 and HCN reach high 
abundances with little dependence on the temperature, once this 
is above a threshold value of about 300 K (necessary for some 
important reactions to overcome their energy barriers). Thus, the 
main differences in the FUV and X-ray models arise because of 
the different dissociating strengths. In fact, if the ionization rate 
is enhanced by a factor of 100 in the X-ray model, then the evo- 
lution of the molecular abundances and their steady state values 
become much closer to those of the FUV model. 

The relative intensity between the FUV and X-ray fields will 
establish whether the chemistry is dominated by FUV photons 
(PDR) or X-rays (XDR). X-rays will drive the chemistry in those 
regions of protoplanetary disks with a column density toward the 
star > 10 22 cm 2 , which a re shielded against FU V photons (A v > 
5) but not against X-rays dMalonev et al.Lfl996l) . In regions with 
column density values > 10 25 cm -2 X-rays are severely attenu- 
ated and do not affect to the chemistry. 

In summary, according to the models the formation of acety- 
lene, hydrogen cyanide and methane occurs in time scales be- 
tween a few years and a few thousands of years, and is clearly 
favored at high temperatures. They reach a steady state in which 
they are destroyed by photodissociation (or X-ray induced re- 
actions) and reactions with H but are continuously re-formed 
through reactions involving H2. Thus, the relative rates between 
re-formation and destruction will determine the steady state 
abundance. A high abundance of H2, maintained due to its con- 
tinuous formation on grain surfaces, is essential for accelerating 
the process 

A+H 2 ^AH + H (1) 

where A is an atom or simple radical. In this way organic species 
AH can reach a high steady state abundance. In fact, if the for- 
mation of molecular hydrogen on grain surfaces is suppressed, 
then H2 is effectively dissociated (due to its continuous partic- 
ipation in reactions of the above type) in about 10 3 -10 4 years 
in the FUV model and 10 6 years in the X-ray model, a time at 
which most of the molecules are also destroyed, thus becoming 
transient species. 

The steady state abundances of C2H2, HCN and CH4 are 
very sensitive to the kinetic temperature (especially in the FUV 
model), and may vary from <10~ 10 at 100 K up to several 10~ 5 
at 1000 K. The chemistry of these simple organic molecules also 
depends, although to a lesser extent, on other parameters of the 
model. The FUV/X-ray field and gas density affect the chemical 
time scale (a decrease in any of these parameters will make the 
chemistry proceed slower) and also to the steady state abundance 
(higher rtu/x or h h/<T ratios favor the formation of molecules 
compared to the destruction by FUV photons or X-rays and thus 
favor a higher steady state abundance). The C/O ratio is also 
an important parameter which affects the wealth of C-bearing 
molecules. Obviously, values close to 1 favor a rich C-based 
chemistry while low values make it more difficult. We neverthe- 
less find that only ad opting C/O ratios su bstantially lower than 
the solar value (0.54; lAsplund et al]|2005l) . then the abundances 
of C-bearing species turn to decrease appreciably. The choice 
of the initial abundances is also an issue when modelling the 
chemistry in protoplanetary d isks and several options have been 
used in the literature (see e.g. Wil lacy et al.ll 998; Sem enov et alj 
120041) . If we assume an initial atomic composition, instead of 
a molecular one, we find that the chemical timescale is notice- 
ably reduced (the step i has been already carried out) although 
the steady state abundance is essentially the same. The non- 
dependence of steady state abundances on the initial composi- 



tion is a consequence of the strong photoprocessing of the ma- 
terial, and indicates that molecular abundances in the PDR of 
inner disks are not affected by the chemical history of the gas. 
The actual situation in PPDs may be different if mixing motions 
bring material from outer cold regions and affect the chemical 
balance in the PDR. 

At this point it would be interesting to discuss our results by 
comparing with some other PDR models appeared in the litera- 
ture. A detailed study of the chemi stry in PDRs has been carried 
out by Sternberg & Dalg arnoldl995l) . They modelled a dense («h 
= 10 6 cm~ 3 ) and highly irradiated (x — 2 x 10 5 ) plane paral- 
lel cloud. They considered a reduced set of chemical species, 
including S- and Si-bearing molecules, of low chemical com- 
plexity (e.g., they did not included species with more than one 
carbon atom). The calculated steady state abundances for C- 
bearing molecules such as HCN and CH4 are lower than 10~ 9 
in the warm region of Ay = 0-2 (7* = 3000-50 K), although rad- 
icals such as CH have abundances as high as 10~ 6 . The most 
remarkable difference between their model and ours lies in the 
rin/x ratio. The much lower value of their model (5 compared 
to the value of 4000 adopted by us) favors photoprocesses com- 
pared to chemical reactions and results in low abundances of 
closed-shell organic mol ecules. A more r ecent PDR model by 
Tevs sier et alj (12004 and lPetv et al] d2005l) . focuses on the study 
of long carbon chains in PDRs such as the Horsehead Nebula. 
They consider a model with a gas density «h = 2 xlO 4 cm -3 
and a FUV field strength x — 60, thus the rtu/x ratio is 300, 
much closer to the value adopted by us. Their model neverthe- 
less predicts C-chains such as C2H and C4H to have abundances 
< 10~ 8 . Besides a different chemical network used, another im- 
portant difference between their models and ours resides in the 
kinetic temperature. In their model the abundance of hydrocar- 
bon radicals peaks in the region of Ay = 1-2, where 7* is lower 
than 100 K. At these relatively low temperatures hydrocarbons 
do not reach high abundances since their formation is inhibited 
by various reactions which have moderate activation barriers. 

4. Abundances in the inner region of a T Tauri disk 

The physical conditions adopted in the chemical models dis- 
cussed in the previous section are assumed to be representa- 
tive of those prevailing in the PDR of inner T Tauri disks. 
Nevertheless, real conditions at different points (r, z) in such 
disks span over a wide range around our adopted values. In or- 
der to check whether photochemistry can account for the forma- 
tion of simple organic molecules with abundances comparable to 
those observed, we have adopted the disk physical structure, gas 
density and temperature at each point, from a st eady state flared 
accretion d isk model provided by P. D' Alessio dD'Alessio et all 
Il998lll999l) . The disk model has a mass accretion rate M = 10~ 8 
M Q yr , grain properties similar to those of interstellar dust 
(maximum dust grain size is 0.25 yum), and a central star with M, 
= 0.7 AT©, r, = 4000 K and R t = 2.6 R . The chemical model 
was run for a grid of 18 x 16 (r x z) points covering the upper 
layers of the very inner disk r = 0. 1-3 AU. The modelled region 
covers the CO/C/C + transition, from an inner height z,„ (where 
the gas is well shielded against stellar and interstellar light and 
all the carbon is as CO) up to the disk surface defined by the 
height Zoo (where the gas is completely exposed to the FUV field 
and all the carbon is as C + , see Fig. [3]). 

We consider that the gas is affected by stellar and inter- 
stellar FUV radiation and by cosmic rays. For simplicity we 
neglect the effect of stellar X-rays, although in the case of 
strong X-ray emitter T Tauri stars the abundances of certain 
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molecules, which are sensitive to the ionization rate such as 
HCO + or HNC, may be substantially affected. The influence of 
an enhanced ioniza tion rate on the molec ular abundances has 
been d iscussed by Aikawa & Herbst ( 1999) and Mar kwick et al.l 
(2002). Concerning the formation of simple organics such as 
C2H2, HCN and CH4, an enhanced ionization rate does not 
greatly change the results from those obtained when only a FUV 
radiation field is considered. 

We treat the FUV field considering that each point (r, z) of 
the disk is exposed to two separate radiation fields, stellar and 
interstellar. The FUV flux emanating from the star is diluted ge- 
ometrically and attenuated by the column density of material in 
the line of sight toward the point (r, z). Interstellar radiation is as- 
sumed to penetrate into the disk in the vertical direction. For sim- 
plicity the radiative transfer of stellar and interstellar FUV pho- 
tons is treated separately, in a plane parallel 1 + 1D approach. The 
computed depth-dependent H2, CO and C photorates are then 
used in the grid of time-dependent chemical models. Although 
a coupling between the attenuation of the FUV field and the 
molecular abundances exists, a self-consistent 2D solution of the 
radiative transfer problem, together with the chemical evolution, 
is beyond the scope of this study. 

The spectral shape of the FUV field emitted by a T Tauri 
star significantly differs from that of the interstellar radiation 
field (ISRF). Therefore, photodissociation/ionization rates can 
be drastically different in the presence of a T Tauri radiation 
field compare d to an enhanced ISRF. This problem has been 
addressed by van Disho eck et al.l d2006l) . who calculated pho- 
todissociation/ionization rates of various species for a 4000 K 
blackbody radiation field, typical of T Tauri stars, and compared 
with the values obtained under the ISRF. They found that, for the 
same integrated flux between 912 and 2050 A, the photorates of 
species such as H2, CO or N2 decrease by several orders of mag- 
nitude when the stellar radiation field, instead of the ISRF, is 
used. 

Given the impact of the true stellar radiation field on the pho- 
tochemistry, we have used the Meudon PDR code to generate a 
new grid of depth-dependent H2 and CO photodissociation rates 
and of C photoionization rates. Instead of adopting a given ISRF 
enhancement, the FUV field is generated by a T Tauri star sim- 
ulated by a blackbody at 4000 K and a stellar radius of 2.6 R G . 
The intensity in the Lyq line, at 12 16 A, can be very important 
in T Tauri stars dBergin et all 120031) . and can have a large im- 
pact on the photorates of some species, such as OH and CH4, 
although it does not affect to species such as H2, CO and N2, be- 
cause their main photoabsorption bands lie short of 1216 A. For 
simplicity, here we do not consider enhanced emission at the 



Lya line. The resulting radiation field at the disk surface is then 
determined by the geometrical dilution from the star. We obtain 
stellar photorates for H2, CO and C that are about 5 orders of 
magnitude smaller than the ISRF photorates (for the same inte- 
grated flux in the FUV), i n agreement with the values reported by 
Ivan Dishoeck et al.l d2006l) . Stellar p hotorates for species other 
than H2, CO and C are taken from Ivan Dis hoeck et al. (2006) 
when available, or are assumed to be equal to the ISRF pho- 
torates otherwise. The visual extinction of stellar and interstellar 
light at a given point is assumed to be proportional to the col- 
umn density of hydrogen nuclei Nh in the direction toward the 
star and outward in the vertical direction respective ly. We use the 
stand ard relation A v = N H (cm 2 ) / 1.87 x 10 21 dBohlin et all 
119781) . adequate for the ISM grain size distribution. In the mod- 
elled region, the stellar FUV field has a strength which is higher 
than the ISRF by orders of magnitude (x = 10 4 -10 7 ). However, 
since the stellar photorates are for many species lower than the 
ISRF photorates by orders of magnitude, the photoprocessing of 
the material in the PDR layers is to a great extent dominated by 
the interstellar, rather than stellar, FUV radiation field. 

The gas densities are 10 8 -10 n cirT 3 in the sampled region 
and the gas kinetic temperature, assumed equal to the dust tem- 
perature in the disk physical model, ranges from 100 to 1000 
K. According to iKamp & Dullemondl d2004f) . the gas tempera- 
ture significantly exceeds the dust temperature in the very upper 
layers with a visual extinction of stellar radiation Ay < 0.1. In 
this superheated surface layer our model most probably under- 
estimates the gas temperature. The implications however for the 
calculated abundances of simple organic species are little since 
this is a region where CO has fully dissociated and most of the 
carbon is as C + . We do not consider adsorption/desorption on 
dust grains since at the high temperatures prevailing in the stud- 
ied layers (7* > 100 K) all the molecules are supposed to have 
been evaporated from grain mantles and to be in the gas phase. 

The abundance distributions at time 10 4 yr, when steady state 
has been reached, for C2H2, HCN, and CH4 are plotted in Fig. [3] 
A thin shell of the disk contains simple organic species with 
abundances ranging from 10~ 8 (in the outer part, at 2-3 AU) up 
to 10~ 4 (in the inner 1 AU where gas is hot). The amount of ma- 
terial contained in such a thin layer is however remarkably large 
due to the large gas densities prevailing there. Vertical column 
densities within the disk PDR can be calculated by integrating 
the gas molecular densities in the z direction (see Fig. @J. The 
column densities of H2 and CO are also shown to visualize the 
molecular abundances relative to these two species. 

lLahuis et alj d2006l) observed toward IRS 46 hot (> 350 K) 
C 2 H 2 , HCN and C0 2 with column densities of 3, 5 and 10 x 10 16 
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crrT 2 respectively. In GV Tau. lGibb et al.l (120071) detected warm 
C2H2 and HCN with column densities of 7 and 4 x 10 16 cm" 2 
respectively, and derived an upper limit fo r methane o f CH4/ CO 
< 0.0037, rather similar to that found by iGibb et alj d2004l) in 
HL Tauri (CH4/CO < 0.005). The inferred warm temperatures 
led these authors to suggest that both acetylene and hydrogen 
cyanide are most likely located in the inner disk, within a few 
AU around the star. In our model, the vertical column densi- 
ties of C2H2 and HCN are a few 10 16 cm" 2 in the inner 1 AU 
but noticeably decrease with r following the decrease in tem- 
perature. The column densities of these two species could be 
enhanced at radii larger than 1 AU if the gas temperatures are 
higher than in our disk model or if vertical and radial mixing 
brings warm material toward the cooler midplane regions, as has 
been suggested to occ ur for CO in the outer regions of PPDs 
( Semenov et all 120061) . Our model predicts CH4 to have a large 
abundance in the very inner disk (< 0.5 AU). However, it is less 
abundant than C2H2 an d HCN throughout the rest of the disk, 
unlike in the model by Ma rkwick et alj (12002) where it is the 
most abundant organic spe cies. The reason s of s uch difference 
are unclear. In the model by Markwick et al. (2002) the large gas 
phase abundance of methane is likely due to its desorption from 
grains surfaces. However, since they do not include formation 
of molecules on grain surfaces, CH4 must have been previously 
formed through a mechanism of gas phase reactions which is 
not explicitly detailed. The regions modelled by them and us are 
however different as we only consider the upper disk layers ex- 
posed to FUV radiation while they consider the full disk in the z 
direction down to the midplane, where most of material is con- 
tained. Thus, the vertical column densities of H2 in our modelled 
region are ~ 10 21 cm" 2 while in their model are ~ 10 25 cm" 2 , so 
that the large abundances calculated by them are coming from 
a region much closer to the midplane than that studied by us, 
which makes difficult a direct comparison between both models. 

In our model, CO2 is formed very efficiently through the re- 
action between CO and OH and has a fairly large column density, 
about 10 1 7 cm" 2 , which agree s with the value observed toward 
IRS 46 bv iLahuis et al.l (l2006). Fig.|4]also gives the column den- 
sities of some other molecules which, although have not been ob- 
served in the inner regions of T Tauri disks, are predicted to have 
relatively large column densities in the PDR. H2CO is formed 
with a moderate abundance through various reactions involving 
atomic oxygen and small hydrocarbon radicals. The presence 
of HNC is strongly related to that of HCN, since it is formed 
through proton transfer to the latter species and subsequent dis- 



sociative recombination. The column density of NH3 has a radial 
profile similar to that of CH4, but with a value about one order 
of magnitude lower. Also it is interesting to note that the abun- 
dance ratio of a closed-shell molecule and its related radical (e.g. 
H 2 0/OH and HCN/CN) is quite large (about 100-1000), in con- 
trast with th e values close to 1 p redicted in the outer disk, r > 
50 AU, (e.g. Will acy et al.l l2000). This is a consequence of the 
large gas densities prevailing in the inner disk which, in spite of 
the strong FUV field, tend to favor the formation of closed-shell 
species through reaction (1). 

5. Conclusions 

In summary, we have shown that a strong FUV/X-ray field may 
efficiently drive a rich C-based chemistry in a dense O-rich gas 
with kinetic temperatures of several hundreds of degrees Kelvin. 
The application of our PDR chemical model to the inner regions 
of a T Tauri disk yields gas phase abundances for simple or- 
ganic molecules such as C2H2 and HCN that range from 5 x 
10~ 5 to 10~ 8 , depending on the radius. This translates into verti- 
cal column densities as large as 10 17 cm" 2 in the very inner disk 
(< 1 AU) down to 10 13 at 3 AU. The model thus explains the 
large column densities of C2H2 and HCN (several 10 16 cm -2 ) 
observed in the inner regions (a few AU) of I RS 46 and GV 
Tau disks dLahuis et al.L 120061: IGibb et alll2007l) . The huge vari- 
ations with radius of the gas density and temperature within the 
PDR layer of protoplanetary disks results in a quite different 
gas phase chemistry, i.e. different molecular abundances and ra- 
tios between them, in the inner terrestrial zone compared to the 
outer regions at several hundreds of AU. Our results await confir- 
mation from further observations of different molecular species 
probing the inner regions of protoplanetary disks. 
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